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ABSTRACT 

Recently, we presented evidence for a cross-correlation of the WMAP satellite map of the 
cosmic microwave background (CMB) and the HEAOl satellite map of the hard X-ray back- 
ground (XRB) with a dimensionless amplitude of 0.14 ±0.05 normalized to the product of 
the rms fluctuations of the CMB and XRB (Boughn & Crittenden, 2004). Such a correlation is 
expected in a universe dominated by a cosmological constant via the integrated Sachs-Wolfe 
(ISW) effect and the level of the correlation observed is consistent with that predicted by the 
currently favored Lambda cold dark matter (ACDM) model of the universe. Since this offers 
independent confirmation of the cosmological model, it is important to verify the origin of the 
correlation. Here we explore in detail some possible foreground sources of the correlation. 
The present evidence all supports an ISW origin. 
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1 INTRODUCTION 



In a remarkably short time, the standard cosmological model has 
changed from a Friedmann universe to a spatially flat universe that 
is dominated by a cosmological constant or some other form of 
dark energy (e.g., Bahcall et al. 1999). So far, the primary evidence 
for this model comes from supernovae redshift/magnitude observa- 
tions (Riess et al. 2004) that imply the expansion of the universe 
is accelerating and from the spatial power spectrum of the fluctua- 
tions in the cosmic microwave background (Bennett et al., 2003). 
While these pieces of evidence are compelling, especially when 
combined with the clustering and density of matter deduced from 
galaxy observations, it is clearly desirable to to seek independent 
confirmation of the ACDM model. 

The late-time integrated Sachs-Wolfe (ISW) effect (Sachs & 
Wolfe 1967) promises to provide such an independent test. In a 
A dominated universe CMB photons undergo a net energy shift 
when traversing linear density perturbations (i.e., 5p/p << 1) at 
relatively low redshifts (z < 1). Crittenden and Turok (1996) sug- 
gested that the ISW effect could be detected by correlating the 
CMB with some nearby (z ;£ 1) tracer of matter, e.g., galaxies or 
AGN. This also occurs in a matter dominated universe with less 
than the critical density; however, the effect is spread over larger 
redshifts (Kamionkowski 1996). Initial attempts to detect such an 
ISW signal (Boughn, Crittenden, & Turok 1998; Kamionkowski & 
Kinkhabwala 1999; Boughn & Crittenden 2002; Boughn, Critten- 
den & Koehrsen 2002) led only to upper limits on A in a flat ACDM 
universe and on the curvature density in an open, matter dominated 
universe. 

The situation changed recently with the release of the first 



year of data from the WMAP satellite (Bennett et al. 2003). Since 
then, six teams have found correlations of the WMAP CMB data 
with tracers of matter including the hard (2—10 keV) X-ray back- 
ground, the NVSS radio galaxy survey, the APM galaxy survey, the 
Sloan Digital Sky Survey, and the 2MASS infrared galaxy survey 
(Boughn& Crittenden 2004a; Nolta, et al. 2004; Myers et al. 2004; 
Fosalba & Gaztanaga 2004; Scranton et al. 2003; Afshordi, Loh, & 
Strauss 2004). All of these results are only marginally statistically 
significant, ;£ 3 o; however, most of them detect an ISW signal with 
an amplitude similar to that predicted by the ACDM model. The 
X-ray/CMB cross-correlation (Boughn & Crittenden 2004) is, per- 
haps, the most statistically significant detection, thanks to the large 
area covered by the X-ray data and the fact the the AGN sources 
are at somewhat higher redshifts. It is the purpose of the present 
paper to provide the details of the X-ray/CMB correlation analy- 
sis and to demonstrate that this correlation is unlikely to be due to 
contamination by foregrounds, such as microwave emission from 
X-ray sources or from Galactic emission. 



2 THE/f£AOlA2 2-10/f£l/X-RAYMAP 

The HEAOl data set we consider was constructed from the out- 
put of two medium energy detectors (MED) with different fields 
of view (3° X 3° and 3° x 1.5°) and two high energy detectors 
(HED3) with these same fields of view (Boldt 1987). The data 
were collected during the six month period beginning on day 322 
of 1977. Counts from the four detectors were combined and binned 
in 24,576 1.3° x 1.3° pixels. The pixelization we use is an equato- 
rial quadrilateralized spherical cube projection on the sky (White 
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and Stemwedel 1992). The combined map has a spectral bandpass 
(quantum efficiency ^ 50%) of approximately 3—17 keV (Jahoda 
& Mushotzky 1989). For consistency with other work, all signals 
are converted to equivalent flux in the 2—10 keV band. 

Because of the ecliptic longitude scan pattern of the HEAO 
satellite, sky coverage and, therefore, photon shot noise are not 
uniform. However, the variance of the cleaned, corrected map, 
2.1x 10~^ (rOJcownfi 5"^ )-^, is significantly larger than the mean 
variance of photon shot noise, 0.8 x 10"-^ {TOT counts 
where 1 TOT counts s^"^ ^ 2.1 x 10^^'erg s^^cm^^ (Allen, Ja- 
hoda & Whitlock 1994). This implies that most of the variance in 
the X-ray map is due to "real" structure. For this reason and to re- 
duce contamination from any systematics that might be correlated 
with the scan pattern, we chose to weight the pixels equally in this 
analysis. However, weighting the pixels in proportional to their sky 
coverage makes little difference in the subsequent analysis. 

The resulting map of the hard X-ray background (XRB) has 
several foreground features that were fit and removed from the 
map: a linear time drift of detector sensitivity, high latitude Galac- 
tic emission, the dipole induced by the earth's motion with respect 
to the XRB, and emission from the plane of the local superclus- 
ter. These corrections are discussed in detail in Boughn (1999) 
and in Boughn, Crittenden, & Koehrsen (2002). The map was 
also aggressively masked by removing all pixels within 20° of the 
Galactic plane and within 30° of the Galactic center. In addition, 
large regions (6.5° x 6.5°) centered on 92 nearby, discrete X-ray 
sources with 2—10 keV fluxes larger than 3 x 10~^'erg s~^cm~^ 
(Piccinotti ct al. 1982) were removed from the maps. Around the 
sixteen brightest of these sources (with fluxes larger than 1 x 
I0~^^erg s~^cm~^) the masked regions were enlarged to 9° x 9°. 

Finally, the map itself was searched for "sources" that ex- 
ceeded the nearby background by a specified amount. Since the 
'quad-cubed' pixelization format lays out the pixels on an approx- 
imately square array, we averaged each pixel with its eight neigh- 
bors and then compared this value with the median value of the 
next nearest sixteen pixels (ignoring pixels within the masked re- 
gions). If the average flux associated with a given pixel exceeded 
the median flux of the background by a prescribed threshold (1.75 
times the mean shot noise), then all 25 pixels (6.5° x 6.5°) were 
removed from further consideration. The result of all these cuts is a 
map with 33% sky coverage. This was the same map that was used 
in previous attempts to detect a correlation of the the XRB and 
CMB (Boughn, Crittenden, & Turok 1998; Boughn, Crittenden, & 
Koehrsen 2002) and is used as our canonical X-ray map in the fol- 
lowing analysis. We have also used other maps with less aggressive 
source removal and the correlation results are nearly independent 
of these cuts; although, the less aggressively masked maps result 
in somewhat higher noise. For details of the various source cutting 
schemes, see Boughn, Crittenden, & Koehrsen (2002). 



3 THE WMAP CMB MAP 

The primary CMB map used in the following analysis is the "in- 
ternal linear combination" (ILC) map derived from the first year 
of data from the WMAP satellite (Bennett et al. 2003) re-pixelized 
to be in the same format as the HEAO X-ray map. The ILC map 
was constructed so as to have little contamination from the Galaxy; 
however, to avoid any residual low Galactic latitude contamination, 
we masked it with the most aggressive WMAP Galaxy mask (kO) 
from Beimett et al. (2003), which results in 68% sky coverage. For 
the 1 .3° angular resolution of the ILC map, instrument noise is neg- 



ligible so we chose to weight each of the unmasked pixels equally 
in all subsequent analyses. To check our results we also used the in- 
dependently derived "cleaned" map of Tegmark, de Oliveira-Costa, 
and Hamilton (2003). It is not surprising that this map yielded re- 
sults that were consistent with the ILC map since the two maps 
were derived from the same primary data. To check for a possi- 
ble frequency dependence of the results we also used Q-band (41 
GHz), V-band (61 GHz), and W-band (94 GHz) WMAP maps with 
the same kO masking. The correlation results of these three maps 
are also indistinguishable from the ILC map as will be discussed 
below. 

4 THE CMB/XRB CROSS-CORRELATION FUNCTION 

A standard measure of the correlation of two data sets is the cross- 
correlation function (CCF), which in this case is defined by 

CCF{Q) = ^l^{Ii-l)iTj-T) (1) 

where the sum is over all pairs of pixels i, j separated by an angle 6, 
/,■ is the X-ray intensity of the ;th pixel, / is the mean intensity, 7) is 
the CMB temperature of the ;th pixel, f is the mean CMB temper- 
ature, and A'e is the number of pairs of pixels separated by 9. The 
CCF(0) generated from the "cleaned" map is consistent with but 
typically ~ 10% larger than that generated from the ILC map. The 
average of these two CCFs is shown in Figure 1. The CCF in this 
figure appears to reveal a considerable level of cross-correlation on 
angular scales 6 ^ 10°. However, this is a bit misleading since the 
error bars in Figure 1 are highly correlated. 

We estimated the standard deviations and correlation matrix of 
the errors in the CCF in two ways. Using the data themselves, we 
computed 400 CCFs by rotating the ILC and "cleaned" maps with 
respect to the X-ray map. By performing two rotations of angles 
larger than 20° , the effects of any intrinsic correlations at small an- 
gular scales, < 10°, are minimized. Of course, there were pixels 
pairs that are repeated in some rotations; however, the number of 
these are very small compared to the number of pixel pairs at each 
angle 9. While the signal distribution in the maps is to a good ap- 
proximation Gaussian, this method enabled us to estimate the errors 
in a way that is relatively independent of the statistical characteris- 
tics of the noise in the maps. As an alternative we generated 1000 
Monte Carlo CMB maps expected from Gaussian fluctuations with 
the ACDM cosmological parameters consistent with the WMAP 
data set. These Monte Carlo maps were cross-correlated with the 
the real X-ray map. We did not generate Monte Carlo X-ray maps 
since the XRB is not as well characterized statistically as the CMB. 
That the standard deviations and noise correlation matrices of these 
two methods agree to within a few percent is an indication that the 
errors are well characterized to this level. Since we did not generate 
Monte Carlo X-ray maps, there was no correlated component in the 
Monte Carlo trial maps and, therefore, cosmic variance of the ISW 
signal was ignored. Because of the large sky coverage, the cosmic 
variance noise of the ISW effect is considerably smaller than other 
noise sources and, in any case, does not change the statistical sig- 
nificance of the detection of a correlated signal in the two maps. 

The CCF due to the ISW effect should be achromatic, i.e., 
independent of the frequency at which the CMB is observed (un- 
like the Sunyaev-Zel'dovich effect). To check this the X-ray map 
was cross-correlated with three WMAP CMB maps at frequencies 
41 ± 4 GHz (Q-map), 61 ± 7 GHz (V-map), and 94 ± lOGHz (W- 
map). These three maps were corrected for Galactic emission using 
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5 THE SIGNIFICANCE OF THE DETECTION 
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Figure 1. The points show the cross correlation observed between the 
WMAP and HEAO-A2 data sets. The heavy solid line represents the pre- 
dictions for the best fit WMAP ACDM cosmology. 



The first three points plotted in Figure 1 (0° ^ 9 ^ 2.8°) are more 
than 2 a greater than zero while the next four points (4.0° ^ 6 ^ 
7.8°) all exceed zero by more than 1 a. However, due to the highly 
correlated nature of the noise, the overall statistical significance of 
the detection is only ~ 3 C as will be discussed below. 

Since each of the values of the CCF is determined from a great 
many pixel pairs (from 7.2 x 10^ / pairs for the 9 = 0° point to 1 .5 x 
10^/ for = 90°), it is reasonable to assume that the noise in the 
CCF is Gaussian distributed (by the central limit theorem). Indeed, 
the observed distributions of the both the 1000 Monte Carlo trials 
and the 400 rotated map trials were consistent with Gaussianity. As 
a test of the tails of the distribution we noted the number of times 
the CCF(O) values of the noise trials exceeded that level observed 
in Figure 1. Out of the 1000 Monte Carlo trials, the noise exceeded 
this level 8 times, i.e., 0.008% of the time. This corresponds to a 
2.4 a effect which is consistent with the value of 2.5 a of Figure 
1. Similar agreement was found for the rotated map trials and for 
other separation angles. 

To determine the overall statistical significance of the detec- 
tion of the ISW effect requires an ISW model. However, the likeli- 
hood that the observed CCF is due to noise alone (with no intrinsic 
correlation of the two maps) can be evaluated in the absence of any 
model, provides such an estimate, i.e.. 



x2 = 52ccF(e,0(Q,r'ccF(e,o 



(2) 
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Figure 2. The cross congelations using the various WMAP frequencies (Q,V 
and W) are all consistent with that found using the ILC map. The nearly 
coincident red, green and blue curves are CCFs of the X-ray map with the 
41, 61, and 94 GHz CMB maps while the points are from Figure 1. Much 
stronger frequency dependence would be expected were the signal arising 
with a typical spectral index expected for radio sources. 



the synchrotron, free-free, and dust maps from the WMAP public 
data set. The resulting CCFs are plotted in Figure 2. If the positive 
correlation reported here were due to radio source contamination, 
one would expect the Q-band CCF to be considerably larger than 
the the W-band CCF, which is clearly not the case. The case against 
radio source contamination is discussed in detail in S6.3 below. 



where (C,j)^' is the inverse of the noise correlation matrix, is 
a conservative estimate of significance since it only assumes Gaus- 
sian noise and nothing about the nature of the signal (other than our 
expectation that the signal occurs at small angles 6 ;£ 10°). In any 
case, depending on the number of data points included in the sum, 
the null hypothesis is excluded at the 93% to 99% C.L., the latter 
of which is for the first 8 data points (9 ^ 9°) with Xv = 19.5/8. 

The ISW effect provides a prediction for the shape and ampli- 
tude of the correlation for a given A. Since the shape varies little as 
A varies, we can use it as a template and fit for the amplitude. Doing 
this, the best fit amplitude of the dimensionless correlation function 
(normalized to the rms fluctuations of the CMB and XRB) at 9 = 0° 
ranges from 0.118±0.049 (2.40 C, 99.2 % C. L., x^ = 1-0/1) for 
a fit to the first two data points to 0. 141 ± 0.048 (2.94 o, 99.8 C.L., 
X^ = 10.9/7) for a fit to the first 8 data points. We take the latter 
as our canonical fit and note that the fit to the first 7 data points 
gives the same result. Fits using from 4 to 8 data points vary by 
only ±6%. 

In the above analysis we have assumed that the noise in both 
the X-ray and CMB maps is Gaussian and this is supported by the 
pixel distribution functions of the X-ray intensity and CMB tem- 
perature. As will be seen below, the distribution of the product of 
signals in pixel pairs is also well fit by that expected for Gaussian 
processes as shown in Figures 3 and 4a. Finally, the frequency with 
which the CCF{Q) generated by Monte Carlo trials exceed a given 
level are consistent with that expected from Gaussian processes out 
to 3 O. The fact that the trials generated from rotated versions of 
the maps agreed with the Monte Carlo trials also indicates that the 
noise is statistically well understood. So from a statistical point of 
view, our 3 O results is well characterized by a 0.999% confidence 
level. 

In the present case, the dominant source of "noise" is the fluc- 
tuations in the CMB from the surface of last scattering and can- 
not be easily eliminated. Even if instrument noise were negligible 
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and the Poisson noise associated witli the discreteness of the X- 
ray sources could be eliminated, one would expect only a factor 
of 2 improvement in signal to noise (Crittenden & Turok 1996). It 
will be some time before the next generation of surveys (e.g., the 
Large Synoptic Survey Telescope and the Square Kilometer Array) 
provides us with new, large-scale mass tracers that will enable this 
improvement to be realized. 



6 POSSIBLE SYSTEMATIC ERRORS 



Since the significance of the detection is limited, it is essential that 
we carefully consider possible systematics which could contribute 
to the cross correlation. 



6.1 Localized contamination 

One important issue is discovering if the signal we observe is due 
to strong correlations between a few contaminated pixels, or is due 
to relatively weak correlations on the whole sky, as predicted by the 
ISW effect. One test of this is to split the data into two hemispheres. 
The fits of the amplitude of the dimensionless correlation function 
to the north and south Galactic hemispheres are 0. 17 ± 0.07 (x^ = 
5.3/7) and 0.11 ±0.07 (y} = 13.8/7) using the first 8 data points. 
While the noise in the detection increases, the signals seen in both 
hemispheres are clearly consistent with each other. 

Another test for localized contamination is to examine the 
distribution function of the individual terms contributing to the 
CCF of Equation 1, i.e., (/,■ — /)(7} — f). To a good approxima- 
tion, the individual pixel distribution functions of both maps are 
Gaussian. If the correlated component of the two maps arises from 
the linear ISW effect, then the values of (/,■ — T) and (Tj — f ) 
should be consistent with a bivariate Gaussian distribution with di- 
mensionless cross-correlation coefficient c(9) = CCF(G) / (Ot-Cx), 
where a\ = {{Ij — I)^) is the variance of the X-ray intensity and 

4 



~ iO^i ^ ^)^) variance of the CMB temperature. (Ox = 

0.124 TOT cnts i and Or = 64juA' for the maps we used.) It 
is straightforward to show (see Appendix A) that the distribution 
function, dN/d^, for the product fj = (/,■ — I)(Tj — f )/aj(5x is 
given by 



dN/dii = 



ir(l-c(9)2)i/2 



(3) 



where Kq is the modified Bessel function, // = /j/(l — c(9)^), and 
A'e is the number of pixel pairs separated by an angle 0. Figure 3 is 
a typical plot of the distribution of pixel pairs (in this case for the 
9 = 2.7° bin), as a function of fj along with the theoretical curve 
given in Equation 3. The error bars are the usual Poisson error bars, 
i.e., where n is the number of pixel pairs in the relevant bin. 
The theoretical curve is not a fit to the distribution function data 
but rather uses the average value of c(9) determined by Equation 
1. The difference between the observed distribution [dN/dfj),,/,,^ 
and the theoretical distribution of Equation 3 {dN /dfj),i,„ for the 
data of Figure 3 is plotted in Figure 4a as fractional residuals, i.e., 
{{dN/dn)obs - {d.N / dij),i,eo) / {dN / dn),heo. For comparison. Fig- 
ure 4b is a plot of the fractional residuals assuming c(9) = 0, i.e., in 
the absence of correlations. The nearly straight line in Figure 4b is 
the theoretical curve with correlations and its slope is a measure of 
the positive cross-correlation, c(9). It is clear from this plot that the 
cross-correlation is manifested for the full range of /j. Were the cor- 
relation due to a limited number of non-Gaussian pixel pairs with 



o 




Figure 3. Each bin of the observed correlation in Figure 1 is the average 
of products of pairs of pixels of the X-ray and CMB maps. This shows the 
distribution of the products for all the pixel pairs contributing to the 2.7° 
bin, normahzed by the rms values of the two maps. This figure and the next 
demonstrate that the observed correlation is not arising from a few pixels 
pairs, but is coming from an asymmetry in the full distribution, precisely as 
expected from a weak congelation of two Gaussian fields. 



large values of fj, one would expect significant deviations from the 
linear behavior exhibited in Figure 4b. 

Equation 1 determined the CCF{Q) from a straight average of 
the correlation arising from the pixel pairs. Alternatively, we can 
fit the distributions (such as Figure 3) for the best dimensionless 
correlation c(9) for each angular bin. To do this, we consider a 
statistic Y.i{{dN /dii)obsiMi) - {dN /d^j),heo{lJi)f /ni where the sum 
is over the the bins at /j/ and n, is the number of pixel pairs in the 
bin, i.e., the Poisson variance. Minimizing this statistic with respect 
to c(9) yields another unbiased estimate for the cross-correlation 
function CCF (9). Figure 5 is the average correlation function of 
such fits to the ILC and "cleaned" CMB maps for a fit range of 
—4.0 < /J < 4.0. The CCF (9) determined in this way is consistent 
with the CCF derived from Equation 1 and depicted in Figure 1. 
The errors in Figure 4 were determined from the same 1000 Monte 
Carlo trials discussed above. Importantly, these estimates were in- 
sensitive to the domains of the fits; i.e., fitting the function over the 
interval — 5.0 < /j < 5.0 was consistent with that for the interval 
— 1.0<//< 1.0. This is evident for the example given in Figure 4b 
where the slope provides a reasonable fit to the data over the entire 
range of /j. 



6.2 Instrumental systematics 

In general, cross-correlations are relatively insensitive to unknown 
systematic instrumental errors, especially if the two data sets are as 
disparate as the HEAO X-ray map and the WMAP CMB map con- 
sidered here. It is difficult to imagine how systematic errors related 
to the instruments could be correlated with each other. To be sure, 
if the systematics are extremely large, then even a small correla- 
tion between them could generate a significant CCF. However, any 
residual systematics in the WMAP data appear to be far below the 
level of inherent fluctuations in the CMB (Bennett et al. 2003). The 
HEAO X-ray map does suffer from a significant linear drift in de- 
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Figure 4. (a) The fractional differences between the observed distribution 
seen in Figure 3 and that predicted for the correlation given in Figure 1 at 
6 = 2.7°. The Poisson error bars are largest at large fj where there are the 
fewest pixels pairs, (b) The fractional differences between the observations 
and the prediction for uncoiTelated Gaussian fields. The heavy line corre- 
sponds to the correlated model used in (a) and the data are clearly consistent 
with the theoretical curve over the whole range of /j's. 



lector sensitivity (Jahoda 1993) but we have fit for and removed this 
effect down to a level below that of the intrinsic fluctuations of the 
XRB (Boughn 1998). Also, the instrument drift results in a signal 
proportional to ecliptic longitude and, therefore, represents a larger 
angular scale structure than that in Figure 1 . Thus we conclude it 
is quite unlikely that instrument related systematics are responsible 
for the observed correlations. 



6.3 Foreground radio sources 

More problematic might be the microwave emission of foreground 
radio sources that also emit in hard X-rays. Certainly radio sources 
are highly correlated with the X-ray background (Boughn 1998) 
and if these sources also emit significantly at microwave frequen- 
cies then a positive correlation between the CMB and the XRB 



Figure 5. The cross-correlation function estimated by fitting Equation 3 
to the distribution function of pixel pairs (e.g., Figure 3) for all angular 
separations. These fits were restricted to pixel pair ampHtudes of ^ 4. 
The results are consistent with the CCF of Figure 1 (soHd line) and imply 
that the correlation is not arising from a few contaminated pixels. 



would be expected. If this were the case, the CCF would have the 
same profile as the X-ray auto-correlation function (ACF). How- 
ever, the CCF profile in Figure 1 is substantially broader than the 
X-ray ACF with a half maximum width of --^ 4° for the former 
and ~ 2° for the latter (Boughn, Crittenden, & Koehrsen 2002). In 
addition, the analysis of the distribution function discussed above 
effectively eliminates the possibility that the CCF signal is due to a 
few "ringers", i.e., a few regions of the sky that are coincidentally 
bright (or dim) in both the X-ray and CMB. 

We can attempt to directly estimate the microwave emission 
of the sources of the hard X-ray background. However, this is dif- 
ficult since recent surveys have not been extensive and are lim- 
ited to relatively low frequencies (< 5 GHz); in addition, radio 
emission is only very roughly proportional to X-ray luminosity. 
Never-the-less, we have made rough estimates by transforming 
the low frequency radio data to WMAP frequencies assuming a 
power spectral index of a --^ —0.4, which is the average observed 
for X-ray selected radio sources (e.g., Reich et al. 2000). Using 
the ROSAT/FIRST analysis of X-ray bright AGN by Brinkmann 
et al. (2000), we obtain a rough estimate of the ratio of the 41 
GHz radio flux density to 2-10 keV X-ray flux for these sources, 
log{Fy/Fx) = —14.2. Combining this with the rms fluctuation of 



the X-ray background, 1.5 x 10 ^ ergs 



(Boughn, 



Crittenden, Koehrsen 2002), the implied correlated component of 
the 41 GHz CMB temperature fluctuations is 8T ~ 0.02 fjK. Such 
a correlated component would result in a dimensionless correlation 
amplitude of 6 x 10^^, a factor of 240 smaller than our observed 
CCF. At the other extreme, from the Chandra analysis of Bauer 
et al. (2002), we deduce a flux ratio of log{F^ / Fx) = —12.0. This 
implies rms temperature fluctuations of 5r ~ 3 /jK and a dimen- 
sionless correlation amplitude of 0.09, only slightly lower than the 
observed. These two estimates bracket those estimated from several 
other observations of the radio/X-ray relation (Ciliegi et al. 2003; 
Georgakakis et al. 2003; Reich et al. 2000). Furthermore, the im- 
plied values at 94 GHz are lower by a factor of (41/94)^ '* ~ 0. 14. 
Even though these estimates indicate that microwave emission from 
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sources is not a dominant problem, the data is not yet good enough 
to make a strong claim. 

A much stronger case against the possibility of radio source 
contamination comes from the achromatic nature of the ISW effect. 
The mean radio spectral index of X-ray selected sources is ~ —0.4 
(Reich et al. 2000) while the blackbody spectral index of the CMB 
is a ~ +2.0 in the Rayleigh- Jeans part of the spectrum. If the CCF 
we observe were due to radio source contamination then one would 
expect the CCF with the 41 GHz WMAP map to be (41 /94)-2-4 ~ 
7 times larger than the CCF with the 94 GHz map. Even inverted 
spectrum sources with spectral indices as large as a ~ +1.2 would 
imply a factor of two difference between these two CCFs. The solid 
and two dashed curves in Figure 2 are CCFs of the X-ray map with 
the 4 1 , 6 1 , and 94 GHz CMB maps while the points are from Figure 
1. It is clear that the difference between them is a few percent at 
most and we therefore conclude that it is extremely unlikely that 
the observed correlation is due to radio source contamination. 

We can also consider how the observed correlations depend on 
the level of point source cuts from the X-ray or CMB maps. Were 
bright point sources a dominant contribution, the cross correlation 
should fall off as the cuts remove more sources. In our canonical 
X-ray map, we avoided potential contamination of nearby point 
sources by masking pixels with excessive X-ray emission (see 
above). The resulting masked map had 33% sky coverage. We also 
considered less stringent cuts leaving more sky coverage (Boughn, 
Crittenden, & Koehrsen 2002), and the observed cross correlations 
were largely insensitive to the level of these cuts. 

As an alternative procedure to remove point sources, we also 
used the ROSAT All-Sky Survey (RASS) Bright Source Catalog 
(Voges et al. 1996) to identify relatively bright sources. While the 
RASS survey has somewhat less than full sky coverage (92%), it 
has a relatively low flux limit that corresponds to a 2 — 10 keV flux 
of ~ 2 X 1 0^ ' erg ' cm^^ for a photon spectral index of a = — 2. 
Every source in the RASS catalog was assigned a 2 — 10 keV flux 
from its B-band flux by assuming a spectral index of— 3<a<— 1 
as deduced from its HR2 hardness ratio. For fainter sources, the 
computed value of a is quite uncertain; if it fell outside the typical 
range of most X-ray sources, — 3 < a < — 1, then a was simply 
forced to be — 1 or —3. It is clear that interpolating RASS flux to the 
2—10 keV band is not accurate, so one must consider the level to 
which sources are masked with due caution. However, we are only 
using these fluxes to mask bright sources and so this procedure is 
unlikely to bias the results. 

We considered maps with ROSAT sources removed at three 
different inferred 2—10 keV flux thresholds. Thirty-four, high 
Galactic latitude RASS sources with fluxes in excess of 3 x 
10^"erg s^^crrT^ were identified and 6.5° x 6.5° regions around 
each source were masked. Recall that this flux level is the nom- 
inal level for the Piccinotti sources that are already masked. The 
resulting map had 51% sky coverage. We also identified and re- 
moved sources with inferred 2—10 keV fluxes in excess of 2 x 
10^ "erg s^^cm^^ (47% sky coverage) and fluxes in excess of 
1 X 10^" erg s^^cm^^ (34% sky coverage). The CCFs of these 
three maps were all consistent with the CCF in Figure 1. Even 
a map with no sources removed other than the bright Piccinotti 
sources (Piccinotti et al. 1982) (56% sky coverage), has a CCF that 
is also consistent with that of Figure 1. Since the observed correla- 
tion is insensitive to the flux level of masked sources, we conclude 
that it is unlikely that the CCF is contaminated by radio emission 
from point X-ray sources. 



6.4 Galactic emission 

It is possible that diffuse microwave/X-ray emission from the 
Galaxy could be a source of correlation; however, the WMAP ILC 
map (Bennett et al. 2003) and "cleaned" map (Tegmark et al 2003) 
were both constructed so as to minimize Galactic emission. In the 
case of the X-ray map, we fit for and removed a small component of 
high Galactic latitude diffuse emission (Boughn 1998). As is true 
for instrument drift, any diffuse, high latitude emission would most 
likely be on larger angular scales than indicated in Figure 1 . Never 
the less, we checked for additional contamination by computing 
the CCF for a variety of Galactic latitude cuts from ±20° to ±45°. 
While the noise in the latter was larger due to lower sky coverage, 
the CCFs for all the cuts were consistent with each other. Finally, 
we can again use the achromatic nature of the observed CCF as ev- 
idence against Galactic contamination. The CCFs of our canonical 
X-ray map with the three WMAP maps (41,61, and 94 GHz) with- 
out corrections for high latitude Galactic emission also agree with 
each other to within a few percent. This would not be the case if the 
CCFs were significantly contaminated by Galactic emission since 
the spectral index of this emission —0.7 > a > —0.3 is so different 
from the spectral index of the CMB, a = 2.0. 



7 INTERPRETATION IN TERMS OF THE ISW EFFECT 

To interpret the CCF of Figure 1 in terms of the ISW effect re- 
quires a cosmological model. We assume a flat, ACDM universe 
with the parameters favored by the WMAP CMB power spectrum, 
i.e., a„ = 0.27, Ha = 0.73,//o = VI km s-^Mpc-^), and « = 1, 
a scale invariant spectrum (Spergel et al. 2003). While this model 
completely determines the ISW effect, the resulting CCF also de- 
pends on the redshift distribution of the X-ray flux, dF/dz, and of 
the X-ray bias, bx, defined by 



where Px is the mean X-ray emissivity, p is the mean density of 
matter and 5 indicates the rms fluctuations of these densities about 
their means. We use the dF/dz of Boughn and Crittenden (2004b) 
that was generated from the X-ray luminosity function of Steffen 
et al. (2003) and Cowie et al. (2003). The bias defined in Equa- 
tion 4 is typically a function of both redshift z and of the scale 
on which the fluctuations are observed. In the present case, how- 
ever, the scales associated with the CCF are large (~ 100 Mpc) 
and the redshifts small (z < 1). In this regime there are good rea- 
sons to believe that the bias is both independent of scale and rel- 
atively small, i.e., bx ^ I (Benson et al. 2000; Fry 1996; Tegmark 
& Peebles 1998). The assumptions of scale and redshift indepen- 
dence combined with dF/dz completely determine the shape of 
the ISW CCF with the amplitude linearly proportional to bx- In 
fact, the shape is relatively insensitive to dF /dz and, therefore, the 
uncertainty in the X-ray luminosity function is relatively unimpor- 
tant. The value of the X-ray bias is taken to be bx = 1.06 ±0.16 
as determined from the clustering of the X-ray background using 
the same ACDM model (Boughn & Crittenden 2004b); this value 
is also consistent with recent observations of QSO and AGN clus- 
tering (Groom et al. 2003; Wake et al. 2004). The theoretical curve 
in Figure 1 is the expected ISW effect for these parameters, i.e., it is 
not a fit to the observed CCF. The %^s of the residuals range from 
0.6/1 for the first datum to 12.4/8 for the first 8 data points. Beyond 
10°, the data are consistent with no correlation. It is clear that the 
observed CCF is consistent with the expected ISW effect. 
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It is possible to perform a maximum likelihood fit of an ISW 
model to the data by allowing one or more parameters to vary, e.g., 
A, dF/dz, bx, etc. Currently, A is quite constrained by the WMAP 
CMB data (Bennett et al. 2003) while, as mentioned above, the 
shape of the ISW signal is relatively insensitive to dP jdz- The X- 
ray bias is, perhaps, the least well known of the model parame- 
ters so we allowed it to vary to minimize x^- The values of bias 
so obtained ranged from fc, = 1.31 ± 0.55 for a fit to the first 3 
data points to = 1.58 ± 0.54 for a fit to the first 8 data points. 
The significance of the detections ranged from 2.4ct to 2.9o, which 
provide another measure of the statistical significance of the detec- 
tion. These values are somewhat higher than but consistent with the 
bx = 1.06 ±0.16 value we derived (Boughn & Crittenden 2004b) 
from the X-ray auto-correlation function. In any case, the errors 
are large and this is certainly not the best method to determine ei- 
ther A or bx- Since the amplitude of the ISW correlation is propor- 
tional to the bias, the above fits provide a useful way of character- 
izing the amplitude of the observed CCF. For example, a fit of the 
bias to the first 8 data points implies a dimensionless correlation of 
CCF /CxCt = 0. 14± 0.05 where Cx and Oj- are the rms fluctuations 
of the XRB and CMB. 

To extract cosmological information requires assumptions 
about the scale and redshift dependence of the bias, as well as 
knowledge of dF/dz- If the bias is constant in both scale and red- 
shift, then the predictions for the dimensionless cross correlation 
will be independent of the bias. However, some uncertainty in the 
theoretical predictions still arise from the uncertainty in dF/dz; the 
amplitude varies about 10-20% for different reasonable assump- 
tions for the redshift distribution. This is to be added to the sta- 
tistical uncertainties of --^ 35%. The fact that the observed value 
is slightly higher than the flyv = 0-73 model predictions suggest 
the data prefer slightly higher A, but not dramatically higher. How- 
ever, these weak detections still allow the exclusion of more radical 
models, where the correlation is expected to be much larger (if the 
matter density is very low) or models where a negative correlation 
is expected, such as very closed models (Nolta et al. 2004). 



8 COMPARISON WITH PREVIOUS RESULTS 

Essentially the same analysis as described above was performed 
in a comparison of the HEAO X-ray map with the COBE satellite 
CMB map (Boughn, Crittenden & Koehrsen 2002). The correlated 
signal evident in Figure 1 is large enough to have been marginally 
detected in that analysis albeit with smaller signal to noise because 
of that map's significant instrument noise and lower angular reso- 
lution. However, we saw no such correlation and proceeded to set 
a 95% C.L. upper limit of the ISW effect at a level lower than the 
positive detection claimed in this paper. To search for the source of 
this discrepancy we constructed a combination map consisting of 
25% of the Q-band WMAP map, 25% of the V-band map, and 50% 
of the W-band map suitibly convolved with the COBE antenna pat- 
tern. This map should be a good approximation to the COBE map 
used in the previous analysis. When the smoothed combination map 
is put through the same "pipeline" as the earlier COBE analysis, 
we detect a statistically significant CCF that is consistent with the 
present result. While the different results for these two analyses is 
a bit surprising it is not totally unexpected. A 2 o noise fluctuation 
in the COBE data would account for the difference. 

Furthermore, it is quite possible that the COBE map contains 
some unknown, low amplitude systematic structure. The CCF of 
the X-ray map with the difference of the COBE and smoothed com- 
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Figure 6. The earlier observations using the four year COBE data showed 
much less correlation, apparently because of noise in the maps. The coiTela- 
tion seen with WMAP (blue), is suppressed when WMAP is smoothed with 
the COBE beam (green); however, it is still much higher than the correla- 
tion with observed with COBE (red). While (COBE - WMAP) should be 
consistent with noise, it is apparently anti-correlated with the X-ray back- 
ground (points with en'or bars). The reasons ai'e uncleai', but could indicate 
a low level systematic in the COBE data. 

bination WMAP maps is shown in Figure 6. It is clear from the 
plot that there is some systematic difference between the two CMB 
maps on angular scales ;£ 10 degrees. Dividing this CCF by the rms 
fluctuations of the X-ray map gives an indication of the correlated 
differences between the two CMB maps, i.e., CCF / Ox ~ 1 to 2 /jA" 
compared to the 71 iaK instrument noise per pixel in the COBE 
map. This level of discrepancy would have gone unnoticed in pre- 
vious comparisons of the two maps. However, it is not possible to 
determine whether the discrepancy is due to such a systematic or is 
simply "unlucky" but statistically possible noise fluctuations in the 
COBE map. In any case, we are led accept the more statistically 
significant result with the much cleaner WMAP data set presented 
here. 



9 CONCLUSIONS 

We conclude that the integrated Sachs-Wolfe effect has been de- 
tected at the ~ 3 O level. There is considerable evidence that the 
detection is not due to spurious systematics or contamination by 
point radio sources. If so, then these and possibly other recent ob- 
servations (Boughn& Crittenden 2004; Nolta, et al. 2004; Myers 
et al. 2004; Fosalba & Gaztanaga 2004; Scranton et al. 2003; Af- 
shordi, Loh,& Strauss 2004 ) offer the first direct glimpse into the 
production of CMB fluctuations and provide important, indepen- 
dent confirmation of the new standard cosmological model: an ac- 
celerating universe, dominated by dark energy. It should be pointed 
out that measurements of the power spectrum of CMB fluctuations 
do not show evidence of increased power on large angular scales 
(6 > 20°) as predicted by the ISW effect, but rather indicate that 
there may be power missing on large angular scales (Spergel et al. 
2003). This deficit is intriguing and may be telling us something 
about the formation of the very largest structures in the universe. 
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The consequences of the ISW effect reported in this letter are pri- 
marily on smaller angular scales and are not in direct conflict with 
the larger angular scale power deficit. 
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APPENDIX A 

Here we briefly derive Equation 3 for the probability distribution 
of the product of two correlated Gaussian distributions. We begin 
by assuming that we have two Gaussian variables, x and y, with 
unit variance and zero mean, and are described by the probability 
distribution: 



P{x,y) 



(27i)det|C|i/2 



exp 



(5) 



The dimensionless correlation is given by the variable c = (xy) , and 
the corresponding covariance matrix is 



1 c 
c 1 



(6) 



The distribution of the product z = xy is given by the marginal- 
ization of the distribution given the constraint: 



P{z = xy) = J dxj dy8{xy-z)P{x,y). 



(7) 



Rewriting the Dirac delta function, 5(xy — z) = l/x5{y ~ z/x) if 
X > 0, we can perform the y integral, leaving: 

' dx. 



P{z) = 2 -Pix.z/x). 
Jo X 



(8) 



Substituting the above form for the two-dimensional distribution 
function, we find 



t(l-c2) 

which can be rewritten as 



ecz/(i-c^) 

31(1 -c2)l 



f°° dx 
/ —exp 
Jo X 



(x2 + \)/(l-c2) 



(9) 



(10) 



If we make the substitutions x = e' and z^ = e^'", the integral can 
be written as 

ecz/(l-c^) 



P{z)- 



3l(l-c2)2 



/■ 



1 



(11) 



Finally, defining t' = 2t — to, the final probability distribution can 
be found: 



P{z) 



31(1 -c2)2 2 



exp[-cosh(f')|z|/(l-c2)] (12) 



■z/(l-c^) 



[ko(\z\/{1-c^)).. 



(13) 



31(1 -c2)- 

where K() is the modified Bessel function. (In evaluating the final 
integral, we have used Eqn. 3.337, Gradshteyn & Ryzhik 1980.) 

This paper has been produced using the Royal Astronomical Soci- 
ety/Blackwell Science M^X style file. 
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